Astronomy & Astrophysics manuscript no. 1564 


February 2, 2008 


(DOI: will be inserted by hand later) 





Jet propagation and the asymmetries of CSS radio sources 

S. Jeyakumar^'^, Paul J. Wiita^, D. J. Saikia^, and Jagbir S. Hooda^ 



o 
o 

> 

o 

(N 
(N 



> 
oo 

in 



o 

Oh; 

6 



^ Physikalisches Institut, Universitat zu Koln, Ziilpicher Str. 77, 50937 Koln, Germany 

Raman Research Institute, C. V. Raman Road, Bangalore 560080, India 
^ Department of Physics & Astronomy, Georgia State University, P.O. Box 4106, Atlanta GA 30302-4106, USA 

National Centre for Radio Astrophysics, TIFR, Post Bag No. 3, Ganeshkhind, Pune 411 007, India 
^ Software Engineering, E*trade Financial Corporation, 4500 Bohannon Drive, Menlo Park, CA 94025, USA 

Received 0000; Accepted 0000 

Abstract. As Compact Steep Spectrum radio sources have been shown to be more asymmetrical than larger 
sources of similar powers, there is a high probability that they interact with an asymmetric medium in the 
central regions of the host elliptical galaxy. We consider a simple analytical model of the propagation of radio 
jets through a reasonable asymmetric environment and show that they can yield the range of arm-length and 
luminosity asymmetries that have been observed. We then generalize this to allow for the effects of orientation, 
and quantify the substantial enhancements of the asymmetries that can be produced in this fashion. We present 
two-dimensional and three-dimensional simulations of jets propagating through multi-phase media and note that 
the results from the simulations are also broadly consistent with the observations. 
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1. Introduction 

Compact Steep-spectrum Sources (CSSs), defined to be 
<20 kpc in size (for i/o=100 km s^^ Mpc'^ and qo = 
0.0) and having a steep high-frequency radio spectrum 
{a > 0.5 where Si, oc have received a great deal of at- 

tention in recent years. The general consensus is that these 
are largely young sources seen at an early stage of their 
evolution (e.g., Carvalho 1985; Fanti et al. 1995; Readhead 
et al. 1996a, b; see O'Dea 1998 for a review). Related to 
the CSSs are the very small (< 1 kpc) Gigahertz-Peaked- 
Spectrum (GPS) sources (e.g., Gopal-Krishna, Patnaik 
& Steppe 1983; Stanghellini et al. 1998). The smallest 
double-lobed radio sources have been christened as the 
Compact Symmetric Objects (CSOs, Wilkinson et al. 
1994; Readhead et al. 1996a; Taylor, Readhead & Pearson 
1996; Owsianik, Conway & Polatidis 1998) which are be- 
lieved to evolve into the Medium-sized Symmetric Objects 
(MSOs, Fanti et al. 1995) and later to the standard large 
double-lobed (FR II) radio sources (e.g. Snellen 2003, and 
references therein) 

High-resolution images of CSSs (e.g. Sanghera et al. 
1995; Dallacasa et al. 1995, 2002 and references therein) 
showed that there is a large range of structures associ- 
ated with CSSs, ranging from double-lobed and triple 
sources to those with complex structures. Saikia et al. 
(1995, 2001, 2002) showed that, as a class, the CSSs 
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were more asymmetric than larger radio sources. They 
also performed statistical tests of the orientation-based 
unified scheme, (Barthel 1989; Urry & Padovani 1995), 
and concluded that there were differences between the 
CSS radio galaxies (RGs) and quasars (QSRs) that were 
consistent with the unified scheme; however, there was 
also substantial evidence for more intrinsic asymmetries 
in CSSs than those which were seen for the larger sources. 
Arshakian & Longair (2000) studied the FR II sources in 
the 3CRR complete sample and also concluded that intrin- 
sic/environmental effects are more important for sources 
with small physical sizes; they also found intrinsic asym- 
metries more important for sources of higher radio lumi- 
nosities. Furthermore, while comparing the hotspot-size 
ratio in radio sources, Jeyakumar & Saikia (2000) found 
that CSSs are generally more asymmetric in this param- 
eter as well; explaining such large ratios of hotspot sizes 
seems to require an intrinsic origin for the size difference. 

Detailed polarization observations of individual 
sources, such as 3C147, show huge differential rotation 
measures between the two oppositely directed lobes, sug- 
gesting their evolution in an asymmetric environment 
(Junor et al. 1999). VLB A observations of 3C43 show 
evidence of high rotation measure where the jet bends 
sharply, suggesting collision of the jet with a dense cloud 
of gas (Cotton et al. 2003). The radio galaxy 3C459 has 
a highly asymmetric structure: the eastern component is 
brighter, much closer to the nucleus and is much less po- 
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larized than the western one, again suggesting interaction 
of the jet with a dense cloud of gas. HST images of the 
optical galaxy shows filamentary structures, suggestive of 
tidal interactions (Thomasson, Saikia & Muxlow 2003). 
Statistically, CSS objects exhibit a higher degree of asym- 
metry in the polarization of the outer lobes at cm wave- 
lengths compared with the larger sources (Saikia & Gupta 
2003; see Saikia et al. 2003 for a summary). These asym- 
metries in the central regions of active galaxies may be 
intimately related to the supply of fuel to the central en- 
gine in these objects, possibly due to interactions with 
companion galaxies. 

All of these studies have indicated substantially greater 
asymmetries for small radio sources. It is clearly worth- 
while to attempt models of their propagation so as to 
see if these differences can be understood. In this pa- 
per we develop a simple analytical model that enables 
us to treat both arm-length asymmetries and luminosity 
asymmetries (Section 2). We then extend the analytical 
model to include relativistic effects so as to allow properly 
for Doppler boosting and time delays if the sources are 
ejected closer to the line-of-sight (Section 3). In Section 4 
we summarize results of our two-dimensional and three- 
dimensional hydrodynamical models of jets propagating 
through an asymmetric dense core, a symmetric galactic 
interstellar medium, and finally, a symmetric uniform in- 
tracluster medium (ICM). We discuss the similar results 
from the analytical and numerical models in Section 5. 



2. Analytical Models for Jets Propagating in an 
Asymmetric Environment 

To estimate the distances out to which jets with equal 
beam powers, but propagating out through different in- 
terstellar media (ISMs), would reach at various times, we 
further generalize the ram pressure balance calculations 
discussed for constant density media by Scheuer (1974) 
and extended to jets leaving a power-law ISM and enter- 
ing intracluster media by Gopal-Krishna & Wiita (1987) 
and Rosen & Wiita (1988). Such models assume that the 
jets start with conical cross-sections and that their heads, 
and associated hot-spots, propagate sub-relativistically. 
We also estimate the ratios of radio luminosities of the 
two lobes as functions of time by expanding upon the cal- 
culations of Eilek & Shore (1989) and Gopal-Krishna & 
Wiita (1991). The latter authors suggested that an "ex- 
cess" ISM in the nuclear regions (central few kpc) , perhaps 
due to the interactions of galaxies, could lead to a larger 
fraction of CSS sources at higher redshifts. Such interac- 
tions could both trigger a phase of jet production from 
the central engine and easily lead to asymmetric distri- 
butions of gas in the inner regions for at least the <10^ 
yr during which the jets propagate through the galactic 
ISM. The propagation through denser media also implies 
higher efficiencies for conversion of jet power to radio lumi- 
nosity (Gopal-Krishna & Wiita 1991; Blundell, Rawlings 
& Willott 1999); this concept is supported by the recent 



findings of Arshakian & Longair (2000) that more lumi- 
nous sources are more asymmetric. 

Our model calculations assume two intrinsically sym- 
metric jets propagating outwards in opposite directions in 
pressure balance between the momentum flux density of 
the beam and ram pressure of the external medium, so the 
equation for the growth of one jet is (e.g., Scheuer 1974), 



dt J ncfimh{D9f 



(1) 



where K is an efficiency factor assumed to be a constant 
and of order unity, finih is the mean particle mass where 
mil is the mass of the hydrogen atom, Lt is the power 
of the beam, D is the distance from the nucleus and 6 is 
the full-opening angle. We assume the gas density of the 
interstellar medium (ISM) to be declining as 



n{D) = r]no 1 + [D/aY 



(2) 



For D ^ a, n{D) « rjno{D/a) This approximation is 
assumed for the analytical model considered here. To allow 
for asymmetrically distributed gas in the central regions, 
the slopes on opposite sides of the nucleus denoted by the 
subscripts 1 and 2, 5i^2 are taken as shallower, and the core 
radius Oc smaller, in the interior regions, but the slope can 
be steeper in the exterior regions. We have considered the 
possibility that the ISM density levels off to a constant 
ICM density at large radii. Standard central densities of 
rio ~ 10~^cm~'^ and a w 2 kpc are typical of large elliptical 
galaxies (Forman, Jones & Tucker, 1985). Excess central 
gas, requiring ry « 10 and extending out to a few kpc 
(i.e., with flc ^ 1 kpc) is sufficient to confine reasonably 
powered jets to CSS dimensions for significant lengths of 
time (e.g., Gopal-Krishna & Wiita 1991; Carvalho 1998). 

We consider a model physically corresponding to jet 
propagation initially through an asymmetric nuclear re- 
gion, followed by a symmetrical galactic scale ISM, and 
finally, through a constant ICM. The central ambient 
medium is characterized by ry = 771 and 5 — 5i on one side, 
with the corresponding values on the opposite side being 
7]2 and 82- For the sake of simplicity, we take the core ra- 
dius for both as ac- The jet propagating outwards on side 
1 reaches a distance Ri after a time Ti, while the one on 
the opposite side reaches R2 after a time T2. Beyond i?i 
and R2, the density distribution on both sides is described 
directly by equation (2), with a 3> dc- The evolution after 
a distance Rz{> Ri) on side 1 and Ri{> R2) on side 2, is 
governed by a constant density, with and r4 the times 
taken by the corresponding jets to reach the distances i?3 
and i?4, respectively (Fig. Although we will usually 
assume R2 = Ri and R4 = R^, the times at which those 
radii are reached will be different for the opposite sides 
even for jets of identical powers, and it is certainly pos- 
sible that R2 ^ Ri or R4 ^ R3, so it is useful to retain 
this more general notation. Quite clearly, the jet propa- 
gating through the denser side (with shallower S or higher 
rj) will take longer to get to any fixed distance. Note that 
there must be a density discontinuity, either at D = if 



S. Jeyakumar et al.: Propagation of jets and the asymmetries of CSS sources 



3 




Fig. 1. Analytic models for jets propagating through asymmetric media: (a) the density distributions, n{D), as a 
function of distance, with side 1 on the right and 2 on the left; since rji — 772 there is a discontinuity at R2; (b) the 
arm-length ratio, ru{t)] (c) the luminosity ratio, rL{t) and (d) the ru — diagram, with both quantities now in log 
units. The open and filled circles represent the times Ti and respectively and the open and filled triangles represent 
the times T2 and T4 respectively. The unit of time is 1 Myr. 



'Hi 7^ ^2, or at i?2 if "Hi = V2] the latter case is shown in 
the left panel of Fig.^ Such discontinuities are most likely 
to arise through an asymmetric distribution of gas from 
a recent merger that may be responsible for triggering 
the nuclear activity (e.g. Gopal-Krishna & Wiita 1991). 
However, it is relevant to note that any realistic situation 
would have steep gradients and not discontinuities. Note 
that even if 5 > 61^2 (as explicitly assumed below), be- 
cause a ^ Uc the absolute rate of decline in density at 
D > i?i,2 is slower than at D < i?i.2 where the inner core 
distributions dominate. 

Using Eq. (2), for D :$> a, equation (1) re- 
duces to D^^-^\dD/dt) = ry-i/V^ifl^--^, where 

i_ 

r = [^nocnmhO^ / 4:K Lb] ^ a? ■ Integrating the above equa- 
tion yields D = [(2 - 5)77-^T-4]i/(2-'5)a. For D = i? at 



t — T, this can be expressed as D = R{t/T) 



1/(2-5) 



We 



use the above expression to estimate the distances Di and 
D2 to the lobes 1 and 2 respectively. 

When the distances of propagation of the jet on both 
sides satisfy Di < Ri and D2 < i?2(= ^^i) then 

Di,2^RiMt/Ti^2f^^'''^'^. (3) 
But once Di > Ri, 

1/(2-5) 

(4) 



2-5 



2-5i 



t 



1 



1 



Ri. 



For i ^ Ti this can be expressed as 

= i?i ([(2 - 5)1(2 - 5^)][tlT,]f'^^''K Similarly, 
the jet propagating on the opposite side, for D2 > R2 

1/(2-5) 

R2, 



D2 = 



F(A-1 
T2 



1 



for 



(5) 



with 



F = 



2-5 



2-52 



52-5 



{aclR2r 



2-5 
2-52 



ac 
Ri 



-5i 



where the last equality holds when rji — 772 and i?i 
For t » T2, D2 = i?2(F/r2)i/(2-^)fi/(2-5). 



R2 



We now estimate the ratios of separations, td, and 
luminosities, r^, for the lobes during different stages in 
the evolution of the source for this model. We choose 
the slopes of the density distributions such that 5 > 
5i > 52, and 772 = rji so that Di is greater than or 
equal to D2 when evaluated at the same time, though 
we leave some dependences on rj explicit in the follow- 
ing expressions. For t < Ti and t < T2, when both 
the jets are within the asymmetric environment, rjj — 
jj.V(2-52)/yi/(2-5i)]^(5i-52)/[(2-5i)(2-52)], Using the ap- 

proximate result that the radio luminosity on each side 
scales as Lr cx n(i:>)3(i+")/i" ex niD)"" (Eilek & Shore 
1989; Gopal-Krishna & Wiita 1991) we find that 

2(5i-52) J 

j,2- 



Hi 

V2 



Ri 



-2S2 

' rp2-S2 

^2 



-4,(61-62) 



where we have taken Ri — R2. With a « 1.0, to « 0.6. 
For the case shown in Fig. 1, where S = 0.75, 5i — 0.675 
and 52 = 0.3, ro oc t°-'^'' , and vl oc t-°-'^°. Thus, at this 
early phase of the evolution, ro increases with time while 
rL decreases with time. 
For t > Ti and t < T2, 



2-5 
2-5i 



m 



2-5 
2-5i 



Ri 



; ^(2-5){2-a2) ; 



2(52— 5i) rp2-5 



-4(S-S2) 
f(2-S)l,2-S2) 



where we have used the continuity of density at Ri to elim- 
inate 77 and a through the relation rya^'' = rjio^.^'^ R^^^^^ 
For the above values of the 5's and m we see that the rise 
in rjj and the decline in rj^ both become slightly stronger, 
with ro oc t°-2i and oc t^°-^^. 

For t > T3, when the faster jet crosses R3, the distance 
travelled by lobe 1 is given by 



61 





( t 








TiJ 



i?3. 



(6) 
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For t T3 this result can be expressed as 

D, = i?3[2/(2 - di)]i{Ri/R,)^-'/'{t/T,)^/\ 

During the last intermediate phase, when t ^ T3 but 
T2 < t < T4, 



and 

rL = 



2-61 



1/2 



R3 

Ri 



^^^2. ^^^,2V(2-^r™ 



To 



For the above choice of the parameters for the density pro- 
file, ru (X i""-^". Thus at moderately late times slowly 
decreases with time, while more rapidly approaches 
unity from below, since tl oc t^ ''^. 

Finally, once t > T^, and the slower progressing lobe 
also enters the uniform intracluster medium, 



Do = 



where 













+ 1 









i?3 



(7) 



02 V'7i 



For t > T4, this simplifies to D2 = i?3G5 (i/ra)^ . Thus, 
when the jets have traversed well beyond R3 and i?4, so 
that t :s> T3,T4, ro — > 1, while ~ 1 as soon as both 
lobes are in the constant density ICM. 

In summary, the size ratio grows monotonically un- 
til T2, after which it declines slowly toward unity. The 
flux-density ratio varies more dramatically, and with more 
complexity, dropping quite a bit until T2 and then rising 
rapidly at that time, because of our assumption of a den- 
sity discontinuity at i?,2- I3oth before and after T^, the flux- 
density ratio of the components becomes symmetric more 
rapidly than does the separation ratio, with oc f^^™, 
but because of its earlier severe deviations from unity, it 
remains more asymmetric. Finally, once lobe 2 also crosses 
the outer interface (at T4) and the constant density ICM 
is entered, the ratio tl very quickly approaches unity, but 
ro does so only asymptotically. 

2.1. Results of basic analytical models 

Fig. n shows the variation of density with distance from 
the nucleus on opposite sides, the variations of the sepa- 
ration ratio, ru, and the flux density ratio, r^, with time 
as the jets propagate outwards, as well as an r/^-r^ di- 
agram. The time axis is in units of 10® yr. The times 
Ti, T2, Ts, and T4 when the jet crosses the interfaces at 
i?2, ^3 J and i?4, respectively, are marked on the tra- 
jectories. Our calculations of such trajectories provide an 
explanation for the observed points in the rjj — dia- 
gram which have ro greater than 1 and smaller than 1 



(see Saikia et al. 1995, 2001). Our calculations also show 
that after the jets have left the asymmetric environment, 
and are passing through similar media on opposite sides, 
the source tends to become symmetric quite rapidly. Thus 
most large FR II sources, which tend to be quite sym- 
metric in both distances (cf. Arshakian & Longair 2000) 
and luminosities of the oppositely directed lobes, may still 
have passed through an asymmetric environment in early 
phases of their evolutions. This is consistent with the pos- 
sibility that most sources may have undergone CSS phases 
during the early stages of their evolutions, as discussed in 
Section 1. In only a few of the CSS sources, such as 3C 
119 (Ren-dong et al. 1991) and 3C 48 (Wilkinson et al. 
1991), does it appear that the jets may have been dis- 
rupted by interaction with the clumpy external medium 
so that they never evolve into large sources (e.g. De Young 
1993; Higgins, O'Brien, & Dunlop 1999; Wang, Wiita, & 
Hooda 2000). 

For a specific reasonable example, consider a jet of a 
beam power of 7 x 10^® erg s~^ and an opening angle 
of 0.04 rad, propagating through a medium with uq = 1 
cm"'^, core radii Oc — 0.45 kpc and a — 2.0 kpc, with 
(5=0.75, (5i=0.675, (52=0.3. The times required to reach 
i?i=5 kpc, i?2=5 kpc, i?3=25 kpc, and i?4=25 kpc are 
0.025, 0.047, 0.19, and 0.22 x 10® yr, respectively The 
above times for a weaker jet of power 2 x lO'*^ erg s~^ , with 
an opening angle of 0.1 radian propagating in a medium 
with smaller Oc = 50 pc and higher central density no = 10 
cm^'^, but otherwise identical parameters, are 0.26, 1.1, 
2.0 and 2.9 x 10® yr, respectively. These ages are smaller 
than any plausible age estimates for larger sources (e.g. 
Leahy, Muxlow & Stephens 1989; Scheuer 1995; Blundell 
& Rawlings 2000). The brief times during which our mod- 
els produce sources which would be classified as GPSs 
(well within i?i ) are consistent with the different analyti- 
cal approximations employed by Carvalho (1998) for jets 
passing through extra dense (but uniform) matter. 

If there is substantially less gas on one side than the 
other, as is the case either if the power-law of the central 
fall-off is substantially less steep on one side than on the 
other, say (5i — (52 ~ 0.4, or if 772 ~ 4ryi, we can obtain fairly 
wide ranges in the predicted ratios of separations and ra- 
dio luminosities, and explain the very asymmetric sources. 
For most sources, a modest difference in S of about 0.15 
turns out to be adequate to produce the observed asym- 
metries, even assuming symmetric jets of equal powers are 
launched from the central engine. 

3. Relativistic and Viewing Angle Effects 

For more detailed comparison with the observations, two 
other important effects must be included: the observed lu- 
minosities are boosted or depressed via possible bulk rel- 
ativistic motions; and the observed positions are affected 
by light travel time effects. 

In this section the subscript 'a' will denote the lobe 
approaching the observer and the subscript 'r' will denote 
the receding component. So Da and D^ are the distances 
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from the core to the approaching and receding lobes, re- 
spectively, while £'1,2 (i) are the distances traveled by the 
jets propagating through the less dense and denser media, 
respectively, as in the previous section; let be the angle 
between the approaching jet and the line of sight to the 
observer. 

In the core's rest frame we are observing at time t. 
The corresponding times at which emission arises from 
the approaching/receding lobe are: ta = t/{\ — (3aC0Sip)\ 
tr = t/ {1 + PrCOSip), where [3 = v/c. Including light travel- 
time effects, where the appropriate velocity is the mean 
speed of advance to that point, the observed distances of 
the approaching and receding lobes are: 



da{t) = A,.(<)sin(^/(1 - 
dr{t) = Dj.{t)smip/ {\ 



[Da{t)/{ct)]cOS^)- 
[Dr{t)/{ct)]cOS^). 



Thus, the distance ratio, including projection and light 
travel time effects is 

da{t) Da{t) l+[Dr{t)/{ct)]cOSip 



(8) 



dr{t) Dr{t) l-[Da{t)/{ct)\cOSif' 

where the superscript * indicates that relativistic effects 
have been included. 

The luminosity corrections also will depend upon the 
Doppler boosting factor, which could be significant if much 
of the emission comes from the hotspots moving at mildly 
relativistic speeds, which is plausible for young sources, 
such as GPS and CSS objects. For a hot-spot or knot 
within a jet, bulk motions produce a change in flux density, 
so that (Scheuer & Readhead 1979) 

Sobsia.r){y) = S,rn{l^)[T{lTf3cOSip)]-^^+"\ 

where the Lorentz factor, F = (1 — /3^)^^/^. The suitable 
values of (3 are the instantaneous rates of advance of the 
lobes. 

In the notation of §2, Sem = La and La (x n{D)™. 
Therefore, the above equations tell us that the observed 
luminosity ratio is: 

Sobs,a{t) _ ^ii.a(ia)[F:(l-/3:COS^)]-(3+") 



rlit) 



Sobs,r {t) 



(tr)[F;;(l + /3*cos^)]-(3+") 



(9) 



_ 1 dDajta) 



p; = i^^^, and F:,, are de- 



Where /^I - 

termined from these instantaneous f3*s evaluated at the 
approaching (or receding) times. To evaluate this pre- 
cisely, albeit in a model dependent fashion, explicit for- 
mulae for Lii,a{ta) and Lii,r{tr) should be used (e.g., 
Eilek & Shore 1989; Blundell, Rawlings & Willott 1999; 
Manolakou & Kirk 2002), and no further simplification is 
possible. However, the boosting term will typically dom- 
inate over the light-travel time factors for the luminosity 
ratio, even for the low values of /3 < 0.2 derived from 
spectral aging and other considerations (Scheuer 1995; 
Gopal-Krishna & Wiita 1996; Arshakian & Longair 2000). 
Therefore, for analytical work Eq. (9) can be approxi- 
mated as 

LR,a{t) [F:(1-/3:cos^)]-(3+-) 



rlit) 



(t) [r;(n-/3,*cos(^)]-(3+") 



(10) 



3.1. Evolutionary tracks in the ro — ri diagram 

Case 1: Here wc consider the situation where the jet 
approaching us is propagating through the less dense 
medium. Then both density and travel-time effects act to 
increase the observed arm- length ratio, rj,. However the 
travel-time and density induced luminosity changes work 
against the relativistic effects in finding the luminosity ra- 
tio, r^. Since the side approaching us comes through a 
lower density medium, it is both farther out and observed 
at a "later" time, so that the density is even lower (or 
at most the same) as it was for the receding lobe; both 
of these effects imply a lower emissivity; however, the ap- 
proaching lobe is Doppler boosted. Meanwhile, the reced- 
ing lobe is still in a higher density region (at least until 
both have reached the ICM), and is observed at an ear- 
lier time in its history, so n{Dr) is even higher, implying 
relatively enhanced emissivity; but it is weakened by the 
Doppler effect. 

Here the appropriate values are Da{t) = Di{t) and 
Dr{t) = D2{t), where Di_2 are obtained from Eqs. (3-7), 
with the appropriate equations used depending upon the 
distances the two lobes have travelled. Then, recalling that 
rD{t) = Di{t) / D2{t) we have for the observed ratio. 



dgjt) 

dr{t) 



rnit) 



l + [D2it)/{ct)]coaip' 
1 - [Di{t)/{ct)]cosip 



(11) 



and 



rlit) ~ rL{t) 



T*2il + l3*2COScp) 

Ft(l-/3*cos^) 



(3+a) 



= rL{t)Vi{t).{12) 



The observed r J, is increased more for smaller viewing 
angles as shown in the upper-left panel of Fig. 2. Although 

can be quite small at very early times, the ratio of the 
Doppler boosting factors, is usually greater than unity, 
and this raises r^. Even for modest /3(< 0.2) it can invert 
the ratio and let the lobe more distant from the nucleus 
appear brighter if it is approaching the observer. At early 
times, the density effects dominate and the sources fall 
below the r*^ — 1 line, unless is very small, in which 
case r£ may initially exceed unity before quickly drop- 
ping below it. Later on, the value of r^ rises above 1 for 
if < 90°, and for most viewing angles the source spends 
a considerable amount of time within a factor of a few of 
the ri^ = \ line, as shown in the upper-middle panel of 
Fig. [3 The evolutionary tracks move to the right on the 
^l ~ ^l) plane as the jet points closer to the line-of-sight 
(Fig. 2, upper-right). Therefore, this scenario can explain 
sources observed to have fairly substantial (~ 10) asym- 
metries in distance and/or luminosity; those asymmetries 
can be correlated or anti-correlated, depending on viewing 
angle and evolutionary phase. 

Case 2: In this situation the side approaching us comes 
through the denser medium, so Da{t) — D2{t) and 
Drit) = Di{t). Now, 



1 l + [Di{t)/{ct)]cos(p 
drit) ^ roit) 1 - [D2it)/{ct)]cosip' 



(13) 
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Fig. 2. Evolutionary tracks for r|j, r£ and — diagrams for different viewing angles to the line of sight. The upper 
panels are for Case 1 while the lower panels are for Case 2. The solid, dashed, dot-dashed and dotted lines indicate a 
viewing angle, (/?, of 10, 45, 75 and 90 degrees, respectively. The symbols marking the time have same meaning as in 
Fig. n and the unit of time is 1 Myr. 
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(14) 



Since l/ri(t) > 1, both of the factors in the above 
expression enhance the luminosity ratio, r^, and at very 
early times it can be very high. However, the two effects 



work oppositely for the arm-length ratio, 



Although 



l/r£)(t) < 1, for smaller viewing angles and at very early 
times, r|j can be well above unity; such sources have 
r'^ > I but are less asymmetric in r|j compared to Case 
1, as shown in Fig. |2| Up until T2, rjj falls; around then 
the density effects are comparable to the light travel time 
effects. At later times r}j rises a bit, but by T4 it starts to 
asymptotically approach unity. Although we have plotted 
r|) (from Eq. 13) as less than 1 in the lower-left panel of 
Fig. 2, since the observed arm-length ratio is defined to 
be > 1, we invert both r|j and for the ip — 90° curve in 
the bottom-right panel of Fig. 2, showing that the evolu- 
tionary track in the r'^ — plane is similar to that of an 
unbeamed source (Fig. 1). A similar procedure is adopted 
for If — 75° in the time range when r|j < 1. The typical 
outcome of this scenario is a source observed to be nearly 
symmetric in ro but possibly very asymmetric in r^. 

4. Numerical Simulations 

We have examined some of our analytical estimates using 
numerical simulations of propagating jets. To follow the 
jets with very high resolution for the extended distances 



that are necessary if we are to examine their evolution 
for situations involving the multiple external media in our 
models, we are constrained to perform two-dimensional 
axisymmetric computations. To allow for some physi- 
cally realistic breaking of axisymmetry, we have also per- 
formed three-dimensional simulations, albeit at lower res- 
olution. We used the well-tested general purpose ZEUS 
2D and 3D (v. 3. 4) codes which are described in depth in 
Stone & Norman (1992a, b). Stone, Mihalas, & Norman 
(1992), Clarke & Norman (1994), and Clarke (1996a). 
As we could not implement the magnetohydrodynamical 
version of ZEUS (e.g., Clarke 1996b), we cannot make 
even the quasi-realistic estimates of the source luminosi- 
ties that can arise from codes that include magnetic fields. 
Therefore, these simulations are useful for examining D{t) 
and ro(t), but not L{t) or rL(t). We now briefly summa- 
rize results from these efforts. 

4.1. 2-D Simulations 

For the first set of simulations we used the ZEUS 2-D 
code as modified by Hooda et al. (1994). Our 2-D work 
used cylindrical coordinates with 900 zones along the axial 
(jet propagation) direction and 200 zones in the radial 
direction. The jet radius has been scaled to 0.05 kpc, with 
the entire grid running from 1-35 kpc in the axial direction 
and 0-2 kpc in the radial direction. With these choices, 
25 zones are included within the initial jet radius, and 
adequate resolution of internal shocks and rarefactions are 
achieved. 




Fig. 3. Two-dimensional simulations, set up to nearly reproduce our analytical model in the plane of the sky {ip — 90°). 
The natural logarithm of the density is shown for Runs A2 (upper image in each panel) and B2 (lower) a,t t — 0.80 
(top panel) and t = 2.05 (bottom panel). The computation is conducted on the upper half of the illustrated region 
and symmetry about the propagation axis is assumed in plotting the density. The total length of the grid is 35 kpc, 
with i?i = 5 kpc and i?3 = 25 kpc; here the unit of time is 0.6 Myr. 



The density within the ISM exactly follows the profile 
as defined by Eqn. (2), and does not make the pure power- 
law assumption used in §2. In this work the important 
parameters are: the jet Mach number, Alj — u/a-m, with 
Oin the jet sound speed and u the velocity of the jet; ^ = 
Pjot/po, with po the initial ambient density and pj the 
density of the jet. Our numerical simulations also assume 
initially conical jets, i.e., we set the half opening angle 
to 0.05 rad in these computations (cf. Wiita, Rosen & 
Norman 1990; Wiita & Norman 1992.) 

We performed two simulations (A2, B2, where A2 cor- 
responds to D2 and B2 to Di of Sect. 2) with parameters 
given in Table 1. Beyond 25 kpc both densities were con- 
stant and equal. In both cases, the jets are assumed to 
be in initial pressure balance with the external medium. 
All individual portions of the external gases are taken as 
isothermal and the interfaces are pressure matched; thus 
in Case A2 there is an abrupt density drop by a factor 
of 0.03 at the first interface at 5 kpc (easily seen in the 
uppermost panel of Fig. , and there is a corresponding 
rise in T by a factor of 33 across that interface. 

Although the code is not relativistic, the simulation 
can be reasonably scaled so that it is applicable to FR II 



Table 1 . The parameters of the numerical runs 





A2 


2D 

B2 


A3 


3D 

B3 


c 


0.001 


0.001 


0.001 


0.001 


Mj 


26 


26 


26 


26 


9j (radian) 


0.05 


0.05 


0.02 


0.02 


i?i(kpc) 


5 


5 


5 


5 


ii3(kpc) 


25 


25 


25 


25 


Be (kpc) 


0.05 


0.05 


0.45 


0.45 


a(kpc) 


1.0 


1.0 


2.0 


2.0 


S 


0.75 


0.75 


0.75 


0.75 


Si 




0.675 




0.675 




0.3 




0.3 




-Ro(kpc) 


0.05 


0.05 


1.0 


1.0 



radio sources. For a value of the initial jet radius i?o = 
0.05 kpc and Tjsm = lO^K, one time unit corresponds 
to 0.6 X lO^yr. The velocity of the jet flow with respect 
to the ambient medium is 0.4c and the power of the jet 
is Lfc ~ 2 X lO''^ erg s^^ for a central ambient density of 
no = 10 cm~^. 
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Fig. 4. The ratios of the location of the Mach disk and 
bow shock of Run B2 to those of A2, r^, as functions of 
time scaled to units of Myr. The triangle denotes the bow 
shock rD and the circle denotes the Mach disk ro. 



Figure 13 shows gas density for the two simulations 
(A2 and B2) at the same early (set of two at the top) 
and late (2.5 times greater; set of two at the bottom) 
times. The early time is chosen to correspond to the point 
when the slower jet (A2) just emerges from the denser 
inner core. At that same time, the identical jet propagat- 
ing through the less dense inner core (B2) has long since 
blasted though its more rapidly declining inner medium 
{Sb2 = <5i > 5a2 = and has gone a good way out 
into the symmetric ISM. So at that time, the ratio of dis- 
tances propagated is nearly maximal, since the asymme- 
try produced by the differences between the central media 
is fully felt. At the late time shown, just before the less 
impeded jet leaves the grid, the ratio of distances is sub- 
stantially reduced, in that the early asymmetry in envi- 
ronment is mitigated by the longer phase of propagation 
through symmetric media. Figure 4 illustrates the 2-D re- 
sults for ru defined using both the Mach disk and the bow 
shock positions; we see that for these simulations there is 
essentially no difference between them, as the jet always 
terminates just a short distance behind the bow shock. 
Since CSS sources are restricted by definition to sizes in 
which galactic halos and any irregularities in them will 
play a large role, the much greater asymmetries observed 
for them (as compared to substantially bigger FR II RGs) 
are easily accommodated. 

4.2. 3-D Simulations 

Three-dimensional simulations allow for the generation of 
non-axisymmetric modes and instabilities. The parame- 
ters of our 3-D simulations are also given in Table 1, and 
are similar to those of the 2-D simulations. The main dif- 
ferences are the need to use a larger values of i?o (1-0 
kpc), Oc, and a and a somewhat smaller value of 0j, 
these changes being necessary to allow sufficient resolution 



in the necessarily coarser 3-D simulations. These three- 
dimensional simulations were run on 150x50x50 active 
zones, extending out to 35i?o (35 kpc) along the initial di- 
rection of propagation, x, and to ±12i?o (±12 kpc) in both 
directions perpendicular to it. To provide marginally ade- 
quate resolution in the most important region, 14 uniform 
zones span the jet diameter, with the rest of the zones log- 
arithmically scaled (e.g. Hooda & Wuta 1996; 1998). This 
set-up allows the jets sufficient time to propagate through 
each medium so that the effects of the different ambient 
media on the jets can be ascertained; see Saikia et al. 
(2003). We believe these to be the first 3-D extragalactic 
jet simulations to employ non-zero opening angles. 

With the exception of a relatively wider jet (in com- 
parison to the size of the core and main galaxy, which is 
required to provide adequate resolution within the jet), 
and minor differences in the scale heights, the parameters 
for the first two 3-D simulations are identical to those of 
the 2-D simulations. Therefore we can quite quickly see 
the differences produced by adding a third dimension. 

Figure El displays densities from simulations A3 and 
B3 at early, middle and late times in the evolutions. For a 
canonical value of i?o = 1 kpc and Tjsm = lO^K, the unit 
of time is now 6.0 x 10^ yr. A major difference between 
the 2-D and 3-D runs is that the bow shocks run ahead of 
the Mach disks by much greater amounts in 3-D than they 
do in 2-D. At early times, when the jet is still in a tran- 
sient initial state, the bow shocks for these Mj = 26 jets 
proceed out through the inner halo very quickly, and open 
a large gap between the bow shocks and the working sur- 
faces at the ends of the jets. Our analytical models assume 
that these two locations stay very close together, so these 
simulations reveal that this assumption is an oversimplifi- 
cation for such high Mach numbers and low jet densities. 
A similar behaviour is also seen in axisymmetric simula- 
tions, using a different 3-D code, of very light cylindrical 
jets propagating through constant density media (Krause 
2003). 

The finite opening angles are fairly well preserved in 
the first phases of the evolution (both in the central halo 
and for some distance into the more steeply declining 
ISM), despite the presence of a recollimation shock. Once 
the jets are at a substantial distance into the ISM, the jet 
becomes better collimatcd and then propagates nearly as 
a cylinder (see the middle panel of Fig. [S] for B3 and the 
right panel for A3). At this point, our analytical models, 
which assume the opening angle stays constant, lose more 
of their validity, as the slowing down of propagation pre- 
dicted by the spreading jet picture is replaced by a situa- 
tion of more nearly constant velocity. We see a slow break- 
ing of axisymmetry through these simulations, despite the 
symmetric initial conditions. This can be attributed pri- 
marily to numerical errors introduced by our setting up 
of a conical ffow within a relatively coarse Cartesian grid. 
Unsurprisingly, the non-axisymmetric effects are stronger 
for A3, which crosses a density jump; such a discontinuity 
helps to induce the growth of perturbations (cf. Hooda & 
Wiita 1996, 1998). 
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Fig. 5. The natural logarithm of the density (with white the densest and black the least dense) is shown for Runs A3 
(top) and B3 (bottom) at three times: t = 0.021, left; t = 0.080, middle; t 0.144, right; here the unit of time is 6.0 
Myr. Note that these are slices through a 3-D computation and no symmetries are imposed, so that non-axisymmetric 
motions are clearly visible in the rightmost panels. The grid runs from 1-35 kpc along the propagation direction and 
±12kpc in the transverse direction. 
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Fig. 6. The arm-length ratios, roit) defined for both the 
bow shock and Mach disk locations of Runs A3 and B3 
as functions of time rescaled to units of Myr. The trian- 
gles denote the bow shock ratio and the circles denote the 
Mach disk ratio. 

Figure El shows the arm-length ratio, r]j{t) for B3 to 
A3; this is plotted using both the bow-shock locations and 
the Mach disk locations. We see that the Mach-disk ratio 
basically follows the trend we would expect, with a sub- 
stantial increase at early times while the faster jet escapes 



the denser core and enters the ISM, and then a more grad- 
ual decline after both jets reach the ISM. Unfortunately, 
our simulations could not be followed long enough to see 
the ratio continue to decline once both jets were in the con- 
stant density ICM, as the faster one leaves the grid before 
the slower one gets to that second interface (Fig.|Bl). The 
bow-shock ratio is less variable, but, after first rising, and 
then slightly declining, it again rises at moderate times 
as the Mach disk in B3 catches up with the bow shock, 
while the Mach disk in A3 is only starting to do so. These 
simulations can be scaled so as to fit FR II radio sources. 
For a temperature of the external gas of about lO^K, the 
velocity of the jet with respect to the ambient medium is 
0.4c and the power of the jet is L{, ~ 7 x 10^^ erg for 
a central ambient density of ng = 1 cm"'^. We have used 
this value of L(, in our analytical calculations (§2 and 3) 
for scaling the time units. 

5. Discussion and Conclusions 

The asymmetric models considered above can explain the 
observed trend in the asymmetry parameters and 
(Saikia et al. 1995, 2003; Arshakian & Longair 2000). Our 
model for jets in an asymmetric environment oriented in 
the plane of the sky explains the large number of galaxies 
and quasars having rjj > I and < 1. Our simulations 
of the jets propagating in asymmetric environments show 
a similar trend in arm-length asymmetry as that produced 
in our analytical calculations. We note that the observed 
ranges oi and can be easily accommodated with the 
kind of asymmetric density profile discussed here. The in- 
clusion of time-delay and Doppler boosting effects can ex- 
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plain the rarer sources with > 1, and the fact that these 
are relatively more common for CSS quasars (presumably 
seen at smaller (p) than for CSS radio galaxies. 

Although our analytical model for radio source evolu- 
tion is certainly oversimplified, somewhat more sophisti- 
cated models (e.g. Kaiser, Dennett-Thorpe & Alexander 
1997; Blundell et al. 1999; Manolakou & Kirk 2002; Barai 
et al. 2004), yield similar trends for distances and radio 
powers with jet power, external medium density and time. 
Larger differences emerge at late times and at high red- 
shifts, when different treatments of adiabatic losses and 
our neglect of inverse Compton losses become important; 
however, these are not relevant for the young CSS sources, 
since the synchrotron losses from these compact lobes will 
always exceed the inverse Compton losses (e.g., Blundell 
et al. 1999). The key advantage our simple model has over 
those more complicated approaches is that we avoid the 
necessity of performing numerical integrals to obtain val- 
ues for Lji, so the estimates of Sects. 2 and 3 are thus 
possible. 

Although the analytical and numerical models agree 
in a qualitative sense, it is not surprising that there arc 
differences between the analytical and the numerical cal- 
culations for the arm-length ratio. Both the 2-D and 3-D 
numerical simulations show a peak value for td oi about 
2.7. At later times the numerically estimated r^'s seem to 
approach asymptotic values greater than 1, although the 
simulations could not be continued long enough to con- 
firm this. The analytical maximum td is about 1.7. The 
time taken to reach the first interface is similar in both 
analytical and numerical calculations, but in the simula- 
tions the second interface is reached later compared to the 
times they are reached in the analytical calculations. This 
implies that even though the calculated and simulated ve- 
locities are similar at early times, the velocity of the nu- 
merical jet, particularly for the 3-D case, is slower than 
our simple ram pressure balance calculations would indi- 
cate. These lower velocities in the numerical simulations 
can be understood because the effective working surface 
in the simulations is somewhat larger than the analytical 
model's jet width, particularly at early times. These slower 
speeds, which are more pronounced for the jet propagating 
through the denser gas, also explain most of the discrep- 
ancies between the maximum values attained by ru, and 
the longer times needed for to decrease towards unity. 

A possible reason why the jet in the 3-D simulation 
slows even more with respect to the analytical estimate is 
that the simulated jet becomes elliptical and rotates about 
its own axis, as expected from theoretical calculations 
of the fastest growing perturbations (e.g. Hardee, Clarke 
& Howell 1995). The jet eventually produces an oblique 
Mach disk, effectively making the jet thrust smaller while 
narrowing the transverse width of the Mach disk (e.g. 
Norman 1996; Hooda & Wiita 1996). Thus, allowing for 
non-axisymmetry in the computations basically acts to 
preserve the asymmetry in m, and could explain the 
larger asymmetry seen in these simulations as well as the 



slower velocities at later times, as compared to the sim- 
plified analytical calculations. 

Recently Carvalho & O'Dea (2002a, b) compared the 
propagation of jets of many different values of ( and 
Mj. These simulations also used the ZEUS 2-D code and 
were restricted to axisymmctry, hut they considered both 
propagation through constant density media and declin- 
ing power laws. Although several parameters are different 
from our study, we made rough estimates of arm-length 
ratios using their high Mj simulations through an ambi- 
ent medium of constant density {S ~ 0.0), and declining 
density {S = 0.75), at a time of 0.5 Myr; td is about 
2.0, 1.3 and 1.1 for C of 0.0276, 0.0542 and 0.1106, re- 
spectively. The ranges of the observed ro thus argue for a 
small C(< 0.01) for extragalactic radio jets. Such low ( val- 
ues have long been known to be favoured by the observed 
widths of cocoons seen in FR H radio sources (Norman et 
al. 1982). 

Carvalho (1998) performed analytical calculations of 
jets "scattering" or penetrating through clouds in an oth- 
erwise uniform medium; the asymmetry is produced due 
to the statistical variance in the number of clouds the 
jets cross through. The arm-length asymmetry produced 
in his model decreases as the jet propagates forward but 
the symmetrization occurs at much smaller scales. To ex- 
plain the corresponding flux density ratios would require 
the shorter arm to be colliding with a cloud at the time 
of observation or at least the presence of more contrast in 
densities between the clouds and regular ISM on the side 
with the shorter arm (Saikia et al. 1995, 2002). In our 
simulations the asymmetry builds up as time progresses 
before getting symmetrized at much later times with the 
corresponding flux density ratio changing in the correct 
sense. 

Recent 3-D simulations of jet/cloud interactions typi- 
cally either show the jets destroying the clouds or being 
halted by them; for a small range of jet/cloud parameters, 
some highly distorted structures and deflected arms can 
be produced (Higgins et al. 1999; Wang et al. 2000). Such 
simulations may be applicable to only a small number of 
CSSs. These cloudy medium models can confine the CSS 
sources to sub-galactic sizes with ages of the order of the 
larger radio sources (Carvalho 1998). But most observa- 
tional evidence (Fanti et al. 1995; Readhead et al. 1996a, b; 
Pihlstrom et al. 2003) suggests that the CSS sources are 
young, not "frustrated" . Our model predictions fit the 
observational studies connecting CSS sources with both 
youth and higher degrees of asymmetry. 

An analysis of the breaks in the spectra of 50 CSS 
sources, which included many measurements at 230 GHz, 
supports the picture that they are very young, with most 
having estimated ages of < lO^yr (Murgia et al. 1999). 
Our models correspond to ages in the observationally es- 
timated ranges of < 10^-lO^yr. Our models also have 
the key advantage of producing the types of asymmetries 
in arm-lengths and luminosities that are characteristic of 
smaller sources. 
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